Partial order in a frustrated Potts model 
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We investigate a 4-state ferromagnetic Potts model with a special type of geometrical frustration 
on a three dimensional diamond lattice by means of Wang-Landau Monte Carlo simulation motivated 
by a peculiar structural phase transition found in /3-pyrochlore oxide KOS2O6. We find that this 
model undergoes unconventional first-order phase transition; half of the spins in the system order 
in a two dimensional hexagonal-sheet-like structure, while the remaining half stay disordered. The 
ordered sheets and the disordered sheets stack one after another. We obtain a fairly large residual 
entropy at T = which originates from the disordered sheets. 
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I. INTRODUCTION 

Generally speaking, frustrated systems have some con- 
straints that forbid simultaneous minimization of all the 
interaction energies. Therefore, frustration usually sup- 
presses phase transitions to long-range orders, and, as a 
result, leads to very rich low temperature phases. More- 
over, a frustrated system may become a spin- liquid phase 
or sometimes exhibit phase transition to a partially or- 
dered state. In this paper, we study a 4-state ferro- 
magnetic Potts model with a special type of geometrical 
frustration on a three-dimensional diamond lattice. This 
model is proposed as a simplified model Hamiltonian for 
a peculiar phase transition found in /3-pyrochlore oxide 
KOs 2 06- Although the obtained phase transition in this 
model does not explain that in KOS2O6, we find that it 
has several interesting properties as a model with geo- 
metrical frustration. 

Let us briefly summarize experimental results of /3- 
pyroc hlore s uperconductors, AOs 2 6 (A is one of K, Rb 
or Cs) 1 * 2 * 3 ^ , from which the particular model studied in 
this paper is derived. We mainly focus on the anomalous 
phase transition other than superconducting transition in 
KOs 2 Og. The pyrochlore structure in general is of type 
A^B^OeO' with cubic space-group Fd3nP^. The four 
crystallographically inequivalent atoms A', B', O, and O' 
occupy the 16c, 16d, 48/, and 86 sites, respectively, in the 
fee unit cell. The crystal structure of the /3-pyrochlore is 
derived from the general pyrochlore structure by remov- 
ing A' atoms, by replacing O' (oxygen) atoms with alkali 
metal atoms, and by filling the B' position with osmium 
atoms. The superconducting transition temperatu re, T c , 
of this /3-pyrochlore oxide family is rather high^^BIIIH 
compared with the previously discovered pyrochlore ox- 
ide superconductor, Cd 2 Re 2 07El. 

Several first-principle density functional calculations 
were perf ormed for the various /3-pyrochlore com- 
pounds^^. The electric structure is not affected by 
changing the alkali atoms which arc almost univalent. 
The osmium atoms and the oxygen atoms form OsOe 
octahedra and OsOe network acts as cages for the alkali 



atoms. The cage size is also insensitive to the alkali atom 
contenlPE Therefore, the difference of T c among the /3- 
pyrochlore family comes from the size mismatch between 
the cages and the alkali atoms. The heaviest anharmonic 
oscillation called "rattli ng" is observed for the smallest 
alkali atom contents, Kj 5 l n l 12 l The spatial asym metry 
of the electron density of K atoms is also observed^ 3 -!, 
which favors the nearest neighbor K atom direction. 

Aside from the superconductivity, the potassium based 
/3-pyrochlorc KOS2O6 undergoes a novel phase transition 
near 7.6KPE3EEI] This is a first-order phase transition 
w here both specific heat and resistivity show hystere- 
g yi6 | i7 | Surprisingly, its transition temperature is in- 
dependent of the superconducting T c . Therefore, this 
transition is not an electric origin and is assumed to 
be of structural one which relates to the "rattling" of 
the K atoms. However, no evidence for lattice dis- 
tortion has been found in X-ray, NMFP^^l high pres- 
sure transport measurements^, scanning tunneling spec- 
troscopy^, the phot o-emission spectroscopy^! or Raman 
spectroscop y^ 3 * 24 ^ measurements. Thus, in spite of the 
very sharp transition, the nature of the phase transition 
remains a big mystery. 

II. MODEL 

Having the peculiar first-order phase transition of 
KOs 2 06 described above in mind, we study a simple clas- 
sical model with special kind of geometrical frustration, 
which was proposed by Kunes et al. using the density 
functional calculation 1 11 1 26 1 27 12 ^. They investigated that 
the on-site potential of the K atoms along with K-K 
bond direction is heavily unharmonic and very flat near 
the symmetry center, which is formed by osmium and 
oxygen surrounding the K ion is over-sized. 

The cages formed by OsOg octahedra in the /3- 
pyrochlore structure have 4 holes towards nearest- 
neighbor K atoms because the 16c sites are empty. There- 
fore these four holes in the cages are not only the origin 
of anisotropic potential of a K atom but also the origin 
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FIG. 1: An approximate length-scale image of the potential 
wells for K atoms. The bond center is the inversion symmetry 
center. 




FIG. 2: (Color Online) Some examples of configurations of 
"bond number", R(ij), and Potts variable a at both ends. 
Each oval means that one of the four wavefunctions, a = 
1, 2, 3, 4, extending in some direction. The green structures 
represent that the Os-O cage have holes towards the next- 
nearest K atoms. Configuration A: two wavefunctions and 
the bond are all parallel. This gives an energy of —2J\ + 2J2 
and is energetically unfavored since J2 > Ji. Configuration 
B: two wavefunctions are parallel but the bond direction is 
different case. This gives an energy of — 2Ji and is energeti- 
cally favored. Configuration C: one of the spins and the bond 
are parallel but the other spins is not. This gives an energy 
of J3. Configuration D gives an energy of 0. 



of rather strong Coulomb interaction between K atoms. 
The K atoms in the /?-pyrochlore KOS2O6 form the di- 
amond lattice structure. Thus the coupling between the 
nearest-neighbor K atoms is essentially of Coulomb ori- 
gin and repulsive. Actually, as shown in the approximate 
length scale of the cage and the flat potential region of 
the K atoms in Fig. [TJ the distance between cages is very 
small. 

For discussing K atoms quantum mechanically, it is 
useful to use the states \a) (a — 1,2,3,4) as the basis 
which represents the sp 3 -like wavefunctions made out of 
Is and 2p states, each located towards the holes of the 
cages. The cages will play the roll of metallic screen- 
ing and thus we assume that the further long-range in- 
teraction can be neglected. The flat potential and the 



holes in cages lessen the energy difference between the 
singlet and triplet states and increases the energy dif- 
ference to the higher energy states. Therefore we focus 
on the lowest four states and neglect the energy differ- 
ence between the lowest singlet and triplet states. In 
this case we can assume the four-fold degenerate basis 
states, I a) (a = 1,2,3, 4), for each K atom in every cage. 
These states can be represented in the sp 3 -like structure 
which point to the nearest-neighbor K atoms. We as- 
sign these states as the four states of the classical Potts 
model. Note that we have one-to-one correspondence be- 
tween the Potts variable 1 to 4 and the direction in the 
diamond lattice. 

We concentrate on the inter-site Coulomb couplings 
represented in the lowest 4 states. Although there 
are 4 4 matrix elements, (af3\W\a' (3') , the largest con- 
tribution comes f rom the diagonal matrix elements, 
{a(3\W\a0 f^' 21 ^ \ These matrix elements arc estimated 
to be 



(a(3\W R \a/3) = -2Ji<J a/9 + 2J 2 6 aK S f } R - J 3 (S aR + 5 0R ) 

(1) 

(Ji = 162K, J 2 = 371K, J 3 = 301K) 



where a and (3 denotes the direction of the bond, 6 a p is 
the Kronecker delta and R is a so-called "bond number" , 
which represents the 4 type of bond directions in the dia- 
mond lattice structure and takes the value of 1 to 4. For 
example, R — 1 means that the direction of the bond is 
(111), and similarly the basis state with a = 1 represents 
one of the four wavefunctions extending in the (111) di- 
rection. 4 types of typical configurations are shown in 
Fig. [2] The first term in the right-hand-side of eq. 
represents the energy gain of 23\ when K atoms on the 
both ends of a bond are in the same direction a = (3. For 
example, the configuration A and B in Fig. [^corresponds 
to this case. The second term in the right-hand-side of 
eq. ([I]) means that there is an energy loss of 2 J 2 when 
the bond direction (R) and the directions of the states \a) 
and |/3) at the both ends of the bonds are all in a straight 
line. The configuration A in Fig. [2] also corresponds to 
this case. The third term in the right-hand-side of eq. ([lj 
represents the energy gain of J3 when the wavefunction 
of one side is parallel to the bond direction. This cor- 
responds to the configuration C in Fig. [2] As a result, 
configuration A in Fig.|2]has an energy, — 2Ji +2J2, con- 
figuration B, — 2Ji, configuration C, — J3 and configura- 
tion D, zero. We can neglect the J 3 terms because these 
terms corresponds to the constant value after summation 
of all bonds. 

The off-diagonal terms, such as (a/3|T4 7 |a7), is of order 
smaller than Ji, J 2 and J3. Therefore we neglect the 
off-diagonalizcd terms in the following. Furthermore, we 
assume that the effects of the excited states higher than 
the lowest four states are negligible. Finally, the effec- 
tive Hamiltonian for the inter-site couplings of K atoms 
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FIG. 3: (Color Online) Illustration of the geometrical frus- 
tration introduced in this model is shown. The first closed 
loop consists of 6 variables and 6 variables ferromagnetically 
orders to "2" state. In this situation, the energy loss in the 
two number "2" bonds, there are no energy loss when we as- 
sign "4" state to this loop but there are direction of the loop 
that "4" state have energy loss. 



becomes 



H = 



(2) 



where R(ij) is also a so-called "bond number" , 1-4, be- 
tween «-site and j-site, which takes the same value of R 
described above. J\ term shows energy gain of parallel 
oscillation of K ion and Ji term shows energy loss of os- 
cillation in line. The bare values of J\ and Ji are 162K 
and 371K, respectively^!. The sum is taken over the K-K 
bond network forming the diamond lattice structure. 

When J 2 = 0, this Hamiltonian describes the classi- 
cal ferromagnetic 4-state Potts model on the diamond 
lattice. The second term gives a peculiar interaction and 
the origin of frustration as shown in Fig. [3} If we consider 
a ferromagnetic state, i.e., a state where all the K atoms 
are in the same wavefunction \a), there must be always 
several bonds where the two wavefunctions at the both 
ends are parallel to the bond direction (Number 2 bond in 
Fig. [3]). This gives a frustration. Note that the term "the 
frustration in the pyrochlore lattice" usually indicates the 
electric one which comes from the tetragonal-structure 
network formed by osmium and oxygen atoms. However, 
the origin of the frustration in KOS2O6 described here is 
completely different. The frustration comes intrinsically 
from the inter-site coupling between K atoms in the dia- 
mond lattice. 



III. METHOD 

The three dimensional 4-state ferro mag netic Potts 
model undergoes the first order transitional In order to 
handle the possible first order phase transition within the 
Monte Carlo simul ation fram ework, we choose the Wang- 
Landau algorithrr J 30 | 31 | 32 | 33 [ This method enables us to 
calculate directly the density of states (DOS), g(E), and 
also allows us to efficiently sample the ground state. This 
algorithm is very effective for studying first-order phase 
transitions. The various thermodynamical quantities are 
also obtained very accurately even near the first-order 
phase transition temperature. Furthermore, this method 



can give estimates of the accurate ground state structure, 
the ground state energy and residual entropy, all of which 
are not obtained before. 

This algorithm works as follows. Since the density of 
states, g(E), is a priori not known at the very beginning 
of the simulation, we first simply set g(E) = 1 for all 
possible energies E. Then we continue to update g(E) 
until it converges to a reasonable functional form and 
the energy histogram h(E) becomes flat. We construct a 
Markov chain of microscopic configuration using the 
local update scheme. We accept the new configuration 
using the transition probability 



POi -> M2) 



(3) 



where E^ and E fll are the energy of the system in a 
specific configuration fii and [i 2 , respectively. The cal- 
culated density of states g(E) and the energy histogram 
h(E) are updated regardless of the acceptance of the new 
configuration as 



In <?(£„) -^IngiE^+hifi 
ft(£7 M ) -> h(E^) + 1, 



(4) 
(5) 



where fi is a modification factor of the ith step of the 
Wang-Landau algorithm, as defined below. At first, /1 
is chosen as f\ = e. After some iteration, we then check 
the "flatness" of the obtained energy histogram h(E) by 
using a criteria that the minimum value of h(E) is not 
less than 80% of the average of the histogram, i.e., 



h{E)\ min > 0.8 x(h{E)) 



(6) 



When the "flatness" of the histogram is reached, go to 
the (i + l)th step of the Wang-Landau algorithm, by 
changing the modification factor as 



ln/i+i = alnfi, 
(0 < a < 1) 



(7) 



where we choose a = 0.5, and also reset the energy his- 
togram. We repeat until the i = 16 Wang-Landau step. 
The choice of i = 16 and a — 0.5 gives a reasonable 
convergence of g(E) for the 4-state Potts modeP^. 

We carry out the Monte Carlo simulations with dif- 
ferent seeds for pseudo random number generator and 
calculated the average and variance of each quantity for 
J2 = 00 and J2 = 10. We use 8 3 = 512 unit cells, each 
consists of 8 atoms and therefore the total number is 
4096. 



IV. RESULTS 

Figure|4]shows the results of energy density for the case 
of J2 = 00. We can see a clear kink at T ~ 0.44 which 
indicates a first-order phase transition of the model. We 
also find that the ground state energy density is equal 
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Temperature 

FIG. 4: Energy density calculation of the extended Potts 
model with J 2 — 00. This results is of 8 3 unit cells and 
consists of 4096 sites. 




Energy Density 

FIG. 5: Weight factor g(E) exp(—f3E) as a function of the 
energy density E, in the case of J 2 = 10 near the transition 
temperature. This results is of 8 3 unit cells and consists of 
4096 sites. 



to —1 per site within the statistical error. If J 2 = 0, 
the ground state energy density should be exactly —2, 
because the model is reduced to the ordinary 4-state fer- 
romagnetic Potts model with 27V bonds in the diamond 
lattice. This result indicates that the ground state of the 
present model is different from the simple ferromagnetic 
ground state due to the frustration induced by J 2 term. 

In order to confirm the first-order phase transition, we 
study the weight factor, g(E) exp(—(3E), as a function 
of the energy density E. Figure [5] shows the obtained 
weight factor of J 2 = 10 and T = 0.435 case. The double- 
peak structure is a clear evidence of the first-order phase 
transition. We confirmed that this double-peak structure 
disappears when the temperature is only slightly off the 
transition temperature, for example, at T = 0.45 ot T — 
0.42. 




FIG. 6: (Color Online) Snapshot of the ground state of 6x6x6 
lattice (1728 atoms). The pole direction of the ellipsoid in- 
dicate the orientation of wavefunction represented by a Potts 
variable.. 

A snapshot of the ground state is shown in Fig. [6] We 
can see that half of the wavefunctions \a) represented 
by Potts variables order in a two-dimensional hexago- 
nal sheet-like structure, in which half sites are located 
slightly above the sheet and the remaining half slightly 
below. This sheet-like plane is perpendicular to the ma- 
jority rattling direction. The sandwiched sheets between 
the directions of the majority wavefunction remain disor- 
dered consisting of the other three wavefunctions. This 
ground state structure is a very new and peculiar ground 
state structure induced by an interesting frustration in- 
teraction in the present model. 

Let us now discuss the residual entropy density origi- 
nating from the disordered sheets. The obtained entropy 
density as a function of T in the case with J 2 = 10 is 
shown in Fig. [7] We find that the entropy density in 
the T — > 00 limit converges to 1.092(6) when we set the 
ground-state entropy to be zero. The correct entropy 
density in the T — > 00 limit should become 

<S)(T^co)=log4, (8) 

since every site has 4 degrees of freedom in the present 
model. Therefore the residual entropy density in the 
present model is calculated as 

(S ) =log4- 1.092 = 0.294. (9) 

This residual entropy density must come from the degen- 
erate ground state of the disordered sheets as shown in 
Fig-i 
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FIG. 7: Entropy density calculations of the extended Potts 
model with J2 = 10 are shown. These results are of 8 3 unit 
cells and consists of 4096 sites and the calculation up to the 
temperature T — 3 is shown. The jump of the Entropy den- 
sity is observed at T ~ 0.44 which corresponds to the first 
order phase transition temperature. The inset shows the cal- 
culation up to the temperature T = 30. The line denotes the 
value of 1.092. 




In the following, we consider the residual entropy in 
detail. Two examples of the ground state configurations 
are shown in Fig. [8] Here the bonds with "bond num- 
ber" 4 are perpendicular to the hexagonal sheet, and it is 
assumed that all the sites in the nearest-neighbor sheet 
have Potts variable 4. Therefore, the Potts variable 4 
is forbidden in the sheets shown in Fig. [HJ Therefore, 
if the Potts variables in the disordered sheets with N/2 
sites are completely random, the total number of possi- 

N 

ble configurations is 3t and the corresponding entropy 
density is 



— log 3": 
N B 



0.5493. 



(10) 



However, the residual entropy density in eq. ^ is much 
smaller than that calculated in eq. ( 10 ). The origin of this 



difference comes from the constraint inside the disordered 
sheet. 

Actually, a bond direction and the two states at the 
both ends of the bond cannot be all parallel inside the 
disordered sheet. Even under the constraint, there are 
many possible configurations in the ground state. The 
upper figure and the lower figure of Fig. [8] look very dif- 
ferent but they have exactly the same energy. The upper 
figure has a site- number unbalance. The ^ sites are filled 
with Potts variable 1 whereas only ^ sites are filled with 
Potts variable 2 and 3. In the lower figure, the sites are 
almost equally filled with the Potts variables 1, 2 and 3. 
Note that the number of the spins in the sheet shown in 
Fig. [8] is not dividable by 3 but there is no difficulty of 
filling Potts variables almost equally. 

The effect of the constraint inside the sheets can be 



FIG. 8: (Color Online) Two examples of the ground state 
configurations in a disordered sheet are shown. Both con- 
figurations satisfy periodic boundary condition. The yellow, 
green and blue ellipse denote Potts variable 1, 2 and 3, re- 
spectively. 



taken into account approximately as follows^ Let us 
start from the completely random configurations which 

N 

consists of 3^ states. If we focus on a single bond in the 
hexagonal sheet, we notice that the forbidden states for 
that bond are included. Since the configurations are cho- 
sen completely randomly, the probability of appearance 
of this forbidden state is |. Here 3x3 = 9 represents 
the total number of the states at the both ends of the 
corresponding bond. Therefore, for each bond, | con- 
figurations should be discarded. As a result, the total 
number of the allowed configurations can be estimated 
as 



3^ 1 



N 



A' 



3JV 



where |iV is the number of bonds in the sheet, 
resulting entropy density is given by 



(11) 



The 



S = log 



= log(1.586) = 0.461, (12) 



which is closer to the numerically obtained residual en- 
tropy density, 0.294, than the completely random value 
0.5493. 



6 



V. DISCUSSION AND CONCLUSIONS 

Let us compare the present results with experiments of 
KOs 2 C>6- Despite the various peculiar properties of the 
present model, it apparently shows directional symmetry 
breaking, which is not observed in KOS2O6. The transi- 
tion temperature of the model, 0.44 Ji ~ 80K, is of order 
high compared with the experimentally obtained first or- 
der transition temperature of T p = 7.5K. Thus, the sim- 
plified model proposed by Kunes et al. does not explain 
the rattling transition in KOs2C>6- A more sophisticated 
model Hamiltonian will be necessary. Recently, Hattori 
and TsunetsugiP^ proposed another model for this rat- 
tling transition of KC^Og- They argued the rattling 
transition by introducing a fifth state in addition to the 
four states discussed above. Although they succeeded to 
explain the first-order phase transition without symme- 
try breaking, the physical origin of the fifth state is not 
clear. Moreover, a recent experiments shows that the 
lattice expands below T pi which contradicts their predic- 
tion. Thus the rattling transition in KOS2O6 remains an 
open question. 

Although the present model is not applied to the tran- 
sition in KOS2O6, we found that this model has very pe- 
culiar features. The obtained phase transition is of first 
order, which is verified from the double peak structure 



of the weight factor shown in Fig. [5] The ground state 
spin configuration of this model shows that half of the 
spins in the system are ordered and form a hexagonal- 
sheet-like structure. The overall possible structures of 
the ground state snapshot configuration are illustrated 
in Fig. [6] and Fig. [8] The remaining half of the spins 
are distributed disordered even in the low temperature 
limit T — > 0, which gives the residual entropy density of 
0.294(6). These results indicate a new kind of interest- 
ing ground state which will be worth further studying as 
a result of new type of frustration. We use the ALPS 
library for calculation^. 
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